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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster  Dimensions 


A  thesis  defense  submitted  in  partial  fulfillment  of  the  requirements 
For  the  degree  of  Master  of  Science  in  Engineering 


Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 


•Objective 

•Section  1:  Introduction 

-  Background 

-  Configuration 

-  Physical  Process 

-  Propellants 

-  Investigation  Constraints  and  Presumptions 

•  Section  2:  Approach  and  Analysis 

-  Discussion  of  References 

-  Analytical  Process 

•Section  3:  Results  and  Discussion 
•Section  4:  Summary  and  Conclusions 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Objective 


The  objective  of  this  study  was  to  develop  the  relationships  required  to  estimate 
small  rocket  thruster  dimensions  over  a  range  of  operating  conditions  as  a 
function  of  the  liquid  propellant  properties  of  viscosity,  surface  tensions,  and 

density. 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Background 

•Alternative  fuel/propellant  sources  being  sought  that  are: 

-  Less  environmentally  impactful 

-  Renewable 

-  Can  be  easily  implemented  into  today’s  infrastructure 

•Fuel/Propellant  candidates  that  meet  these  criteria  will  most  likely  have  different 
formulations  than  those  that  are  the  state  of  the  art  in  common  use 

-  Formulation  changes  can  drive  changes  in  energy  conversion  devices  in  common 
use 

•  Estimation  of  the  delivered  performance  prior  to  empirical  development 
reduces  the  cost  implementation 

•  The  instant  study  considers  how  the  changes  of  a  liquid  propellant’s  physical  properties 
can  impact  the  dimensions  of  an  energy  conversion  device 

-  In  this  case  study,  for  a  spacecraft  thruster 

-  Hydrazine,  N2H4,  and  hydrogen  peroxide,  H202,  are  contrasted 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Configuration 


•  A  two-stage  rocket  thruster 
configuration  was  used  in  the  study 

-  5  lbf  (22  N)  thrust  parameter  was 
set 

•  Used  to  set  the  overall  mass 
flow  rate 

-  Propellant  is  injected  in  two 
locations 

•  First  stage-some  fraction  of 
propellant  is  injected  into  an 
ignition/ decomposition 
device 


Second  ary  Propellant 
Inject  ion  Points 


Energy  conversion  from 
liquid  propellant  to  hot 

_ 


Liquid  Propellant 


Liquid 

Propellant 


Liquid  Propel Lant  Sheet  resulting 
from  impingement  of  liquid  jets, 
disintegrating  into  droplets 


Liquid 

Propellant 


•  Second  stage-remainder  of  propellant  is  injected  downstream  into  the  high 
temperature,  low  velocity  gas  flow 

-  Two  impinging  streams  that  impact  and  atomize  into  droplets 

-  The  high  temperature  environment  prompts  the  droplets  to  convert  from 
liquid  phase  to  decomposed  gas  phase 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Physical  Process 


•  Two  stage  thruster  configuration  dimensions  are 
driven  by  multiple  physical  processes 

-  Stage  1 :  Propellant  decomposition  to  form  a 
high  temperature,  low  velocity  environment 

*  Propellant  decomposition/ignition  means  is 
presumed  to  be  a  catalytic  reactor 

-  Technology  is  established  and  flight 
proven 

•  The  decomposition  species  and  gas 
temperature  exiting  and  propellant  mass 
flow  rate  into  the  decomposition/ignition 
means  are  tracked 


Energy  conversion  from 
liquid  propellant  to  hot 


Liquid  Propellant  I  % 


High  temperature,  low 
velocity,  flow  field 


Stage  1:  Propellant  Decomposition/Ignition 


—  Decomposition  species  knowledge  is  relevant  to  determining  environmental  parameters 
»  specific  heat  ratio  &  heat  transfer  coefficients 

-  Gas  temperature  is  derived  from  specific  impulse  calculations  and  inputted  into  relations 

-  Propellant  mass  flow  rate  into  the  decomposition/ignition  means  drives  atomization 
properties  of  the  second  stage 

•  Propellant  reaction  kinetics  and  species  properties  are  primary  configuration  drivers  from  this 
stage 

-  Establishes  the  environment  that  liquid  propellant  is  injected  into  and  atomized  within 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Physical  Process 


Two  stage  thruster  configuration  is  driven  by 
multiple  physical  processes 
-  Stage  2:  Downstream  Propellant  injection 

•  Propellant  is  injected  as  two  impinging 
streams 

-  Impinging  streams  form  a  liquid 
sheet  which  then  deteriorates 
into  ligaments,  then  droplets 

•  Propellant  physical  properties  of 
density,  viscosity,  and  surface  tension; 
configuration  parameters  of  injector 
orifice  diameter  and  injection  velocity 
are  the  primary  chamber  length 
impacts  in  the  relations  governing  the 
second  stage 


Liquid 

Propellant 


Stage  2:  Propellant  Atomization 


Liquid  Sheet  Deterioration  to  Droplets 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Physical  Process 


Low  jet  velocity  (6.4  m/s)  turbulent  impingement 
sheet  photo  (figure  3a  from  Ryan) 


Low  jet  velocity  (7.1  m/s)  laminar  impingement 
sheet  photo  (figure  4a  from  Ryan) 


•  Photos  of  turbulent  and  laminar  liquid  sheet  deteriorations  into  droplets 

•  Images  taken  from  the  primary  reference,  Ryan  et  al.1 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Propellants 


•  Hydrazine,  N2H4,  and  Hydrogen  Peroxide,  H202  were  propellants  considered 

-  Represent  current  state  of  the  art  propellants  used  for  spacecraft  chemical  propulsion 
applications 

-  Significantly  different  from  each  other  in  physical  properties  and  energy  content 


Propellant 

Liquid  Property 

Symbol 

Value  Range 

Units 

Value  Const!  aints 

Coefficient  of 
Viscosity , 
Dynamic 

Pl 

9.736 X  Iff4 

kg/(m-s) 

2(fC 

Coefficient  of 
Viscosity , 
Kinematic 

9.654X10  7 

m/s 

2(fC 

Density 

Pl 

L009XW 

kg/m1 

2(fC 

Surface  Tension 

°>. 

6.850  X 10" 

kg/s 2 

2<fC 

h2o2 

Coefficient  of 
Viscosity , 
Dynamic 

Pl 

1.250X103 

kg/(m-s) 

2(fC 

Coefficient  of 
Viscosity , 
Kinematic 

VL 

8.621  X 107 

m/s 

2(fC 

Density 

A 

1.450  X  Iff1 

kg/m J 

2<fC 

Surface  Tension 

°L 

8.013X10^ 

kg/s 2 

2<fC 

Summary  of  propellant  physical  properties  of  hydrazine  and  hydrogen 

peroxide 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 


Section  1:  Introduction/Propellants 


•  Reaction  Kinetics 
-  Hydrazine,  N2H4 

•  3  N2H4  4  (1-X)NH3  +  (1+2X)N2  +6XH2 

-  Where  X  represents  the  amount  of  ammonia,  NH3,  dissociation 

•  At  50%  ammonia  dissociation 

-  3  N2H4  — >  2NH3  +  2N2  +  3H2 

-  TC=1200-1300K 

•  At  100%  ammonia  dissociation 


3  N2H4 


3N2  +  6H2 


-  TC=874K 

*  Increased  ammonia  dissociation  leads  to  decreased  gas  decomposition  temperature 
(decreased  chamber  temperature) 

-  Hydrogen  Peroxide,  H202 


•  2  H202 


2  H20  +  02 


-  TC=1274K 

-  Presumes  100%  concentration 

-  Reduced  concentration  leads  to  decreased  gas  decomposition  temperature  (decreased 
chamber  temperature) 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Investigation  Constraints  and  Presumptions 

•  Weber  (We)  and  Reynolds  (Re)  Number 

-  The  relationships  from  Ryan  et  al.1  to  determine  liquid  impingement  sheet  and  atomized  drop 
dimensions  were  bounded  by  Weber  and  Reynolds  Number  ranges  of  validity 

•  350<We<6,600 

•  2,800<Re<26,000 

•  Specific  Impulse 

-  Specific  impulse  values  were  evaluated  across  a  range  of  operating  conditions 

•  The  vacuum  operating  condition  with  an  expansion  ratio  of  50:1  was  chosen 

-  Showed  the  least  variation  across  the  chamber  pressure  region  of  interest 

•  Discharge  Coefficient 

-  Median  value  of  0.77  was  used  in  calculations 

•  Range  of  0.732  to  0.809  resulted  from  calculations  across  the  valid  We  and  Re  range 

•  Injector  Orientation 

-  Impingement  half  angle,  0,  of  40°,  was  chosen  for  calculations 

•  Impingement  Sheet  Dimensions 

-  Angular  displacement,  (|),  from  the  impingement  point  on  the  liquid  sheet  was  presumed  to  be  0° 

•  The  item  of  interest  where  this  is  relevant  was  the  impingement  sheet  length,  which  would 
be  at  (|)  =  0° 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  1:  Introduction/Investigation  Constraints  and  Presumptions 


•  Droplet  Formation,  Size  Distribution,  and  Decomposition 

-  The  effects  of  secondary  atomization  were  not  considered 

•  Eases  computational  complexity  of  the  model 

-  The  droplets  evolved  from  the  atomization  that  will  have  the  most  impact  on  chamber  length  are 
those  evolved  from  the  very  end  of  the  impingement  sheet 

•  Anticipate  these  drops  will  be  the  largest  since  at  the  longitudinal  end  of  the  sheet  the  drop 
velocity  will  be  the  lowest 

-  Lower  sheet  velocities  result  in  larger  droplet  diameters 

-  Decomposition  process  was  considered  to  have  been  achieved  at  the  point  in  time  when  the 
propellant  droplet  had  transitioned  form  a  liquid  to  vapor  state 

•  Vapor  Phase  Thermal  Conductivity 

-  Vapor  phase  thermal  conductivities  were  calculated  by  combining  the  sum  of  reaction  product 
species  by  mol  fraction  presence 

-  More  complex  methods  exist  to  derive  this  quantity,  set  forth  by  Saxena6  and  Mason7 

-  Hydrazine  vapor  phase  conductivity  data  is  not  readily  available,  used  a  surrogate  molecule- 
Methanol,  CH3OH 

•  Hydrazine  vapor  is  not  stable  in  the  temperature  regimes  the  value  is  required  from 

•  Hydrogen  peroxide  data  was  attainable 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Discussion  of  References 


•  Ryan  et  al.1 

—  Empirical  study  that  resulted  in  algorithms  for  estimating  impingement  sheet  length  and  evolved 
droplet  diameter 

•  Lefebvre5 

-  Injector  orifice  estimation  algorithms  relating  liquid  properties  to  discharge  coefficient 

•  Turns9 

—  Evolved  droplet  transition  time  relations  for  liquid  to  vapor  phase  that  consider  environmental 
conditions 

•  Sutton11  and  Hill  et  al.12 

-  Standard  equations  relating  various  aspects  of  rocket  engine  performance 

•  Schmidt2,  Brown3,  The  Hydrogen  Peroxide  Handbook 4,  Svehla8,  and  The  CRC  Handbook  of  Chemistry 
and  Physics 10 

—  Physical  property  information  for  propellants 

•  Spalding13,  Anderson14,  Ibrahim15,  and  Dombroski16 

-  Relevant  citations  from  Ryan  and  Hill  in  regards  to  combustion  chamber  dimension  estimation 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 


•  Operational  conditions  considered 

-  3  Combustion  chamber  pressures  considered 

•  125, 250,  and  500  psi 

-  3  mass  flow  rate  ratios  considered 

•  m’gg  is  the  fraction  of  propellant  allocated  to  the 
decomposition/ignition  device,  stage  1 

•  m’inj  is  the  fraction  of  propellant  allocated  to 
injection  and  atomization  downstream,  stage  2 

•  m’gg  to  m’inj  ratios  are  1:3, 1:4,  and  1:5 

-  2  pressure  “hardness”  ratios  considered 

•  Ratio  of  liquid  propellant  feed  pressure  to 
chamber  pressure 

•  1.25  and  1.75 

•  Indexing  used  to  track  each  case  through  calculations 

-  Combustion  chamber  pressure.Mass  flow  rate 
ratio.Hardness  ratio 

-  Shown  in  table 


Index 

Pc  (psi) 

m'  Ratio 

hardness 

ratio 

l.Li 

125 

0.33 

1.25 

1.2.1 

125 

0.25 

1.25 

1.3.1 

125 

0.20 

1.25 

1.1.2 

125 

0.33 

1.75 

1.2.2 

125 

0.25 

1.75 

1.3.2 

125 

0.20 

1.75 

2.1.1 

250 

0.33 

1.25 

2.2.1 

250 

0.25 

1.25 

2.3.1 

250 

0.20 

1.25 

2.1.2 

250 

0.33 

1.75 

2.2.2 

250 

0.25 

1.75 

2.3.2 

250 

0.20 

1.75 

3.1.1 

500 

0.33 

1.25 

3.2.1 

500 

0.25 

1.25 

3.3.1 

500 

0.20 

1.25 

3.1.2 

500 

0.33 

1.75 

3.2.2 

500 

0.25 

1.75 

3.3.2 

500 

0.20 

1.75 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm- Valid  Cases 

•  Injector  discharge  coefficients  were  first  calculated  using  relations  from  Lefebvre5  based  on  the 
range  of  validity  for  Re  from  Ryan  et  al.1 


CDmax=0-827-0.0085 


k 

do 


1 

Q> 


1  20  1* 
c^V*2'25^ 


•  Total  mass  flow  rate  (m’total)  derived  from  thrust  (FT)  and  specific  impulse  (Isp) 


-  _  Ft 

“““'IspS, 

•  Mass  flow  rate  of  propellants  through  the  injectors  (m’inj)  derived  from  flow  rate  ratios 

-  nigg 

mtotai(l  -  — “  ) 
minj 

m  = - 

inj  #  of  injectors 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm- Valid  Cases 


•  Injector  orifice  dimensions  (A0,  d0)  derived  from  injector  mass  flow  rates,  discharge  coefficient  (CD), 
propellant  density  (pL),  and  pressure  drop  (AP) 


A.= 


Cn\/2pAP 


•  We  and  Re  calculated  from  injection  velocity  (Uinj,  derived  from  m’inj),  orifice  diameter,  propellant 
density,  and  propellant  viscosity  (|lil)  and  surface  tension  (aL) 

-  Cross  check  ranges  of  validity,  per  Ryan  et  al.1 


U.  ^ 

rn 


PtA° 


w, 


PlU 


2  d 

fflj  o 


Or 


Pj 

h 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm-Chamber  Environment 


•  Calculate  specific  heat  ratio  (y)  of  the  reaction  products  as  a  function  theoretical  decomposition/combustion 
temperatures  from  tabular  data  and  ideal  gas  relations  on  mol  fraction  basis 


Ru~ Cp-cv 


Y 


Hi 


vr 


Tirr 


Ym«C-L> 


£y; 


cp(T) 
R 


u 


<*(T> 


•  Calculate  chamber  gas  density  (pc)  based  on  chamber  pressure  (Pc)  and  Temperature  (Tc)  using  ideal  gas 
relations 


Pc= 


Pc 


R«pecTc 


•  Calculate  chamber  gas  to  propellant  liquid  density  ratio  (S) 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 


Section  2:  Approach  and  Analysis/Analytical  Process 


General  Algorithm-Impingement  Sheet  and  Drop  Dimensions 


•  Calculate  sheet  velocity  (Us) 

Us=UjjjjCos(0) 

•  Calculate  wave  number  (km)  from  chamber  gas  density,  impingement  sheet  velocity,  and  surface 
tension 


•  Calculate  impingement  sheet  length  (rb)  from  orifice  diameter,  chamber  gas  to  liquid  propellant 
density  ratio,  and  Weber  number 


-2  2. 
2(14.2S3Wj) 


•  Calculate  sheet  thickness  (h)  from  orifice  diameter,  impingement  sheet  half  angle  (0),  sheet  angular 
displacement  (<|)),  and  impingement  sheet  length 


h= 


d^sin30 


4rb(l-cos<|>cos0)2 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm-Impingement  Sheet  and  Drop  Dimensions 


•  Calculate  liquid  sheet  maximum  growth  rate  factor  (ft  m)  from  kinematic  viscosity  (vL),  wave  number, 
chamber  gas  density,  sheet  velocity,  surface  tensions,  propellant  liquid  density,  and  sheet  thickness 


Pi>m= 


Ulkn 


(-1  + 


1  + 


8(pckmUl-qLk^ 

^LkmPLh 


•  Calculate  drop  diameter  (dD)  from  orifice  diameter,  sheet  velocity,  Weber  number,  chamber  gas  to 
liquid  propellant  density  ratio,  sheet  length,  and  liquid  sheet  maximum  growth  rate  factor 


2.62 

& 

U/b 


1 


J3S6[wcf(0)j3 


m- 


(1-COS0)2 
sin3  ft 


i 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm-Droplet  Evaporation  Quantities 

•  Calculate  chamber  environment  thermal  conductivity  (kinf)  from  tabular  data  and  mol  fraction  for 
reaction  products 


kmf^y.k.OTbar) 


•  Calculate  median  chamber  temperature  (Tbar)  from  propellant  boiling  point  (Tboil)  and  chamber 
temperature 


_(Tc+Tboil) 
bar  2 

•  Calculate  overall  chamber  gas  phase  thermal  conductivity  (kg)  from  propellant  vapor  phase  thermal 
conductivity  (kF)  and  chamber  environmental  thermal  conductivity 

kg=0.4kF(@  Tbar)+  0.6kM(@  Tbar) 


. . - 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 


General  Algorithm-Droplet  Evaporation  Quantities 


•  Calculate  Spalding  number  (Bq)  from  specific  heat  at  constant  pressure  for  gas  phase  (cpg),  chamber 
temperature,  propellant  boil  point,  and  latent  heat  of  vaporization  (hfg) 

-  ratio  of  heat  capacity  as  a  function  of  the  temperature  difference  between  the  environment  and 
propellant  boil  point  to  the  latent  heat  of  vaporization 

T»  _Cpg(TC"Tboil) 

N  x. 

hfR 

•  Calculate  overall  evaporation  constant  (K)  from  the  overall  chamber  gas  phase  thermal  conductivity, 
Spalding  number,  propellant  liquid  density,  and  specific  heat  at  constant  pressure  for  gas  phase 

K=  ln(Bq+l) 

Pi  Cpg 


' 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  2:  Approach  and  Analysis/Analytical  Process 

General  Algorithm-Chamber  and  Throat  Diameter 


•  Calculate  sonic  throat  area  (A*)  from  specific  heat  ratio,  propellant  mass  flow  rate,  chamber 
pressure  and  temperature 


m 


T 

cc  *  c 


7+1 

2  2JFI5 


•  Calculate  chamber  cross  sectional  area  (Achamber)  from  sonic  throat  area,  chamber  mach  number  (M, 
presumed  to  be  0.1),  and  specific  heat  ratio 


Chamber 


—  2  ( 


Calculate  chamber  diameter  (Dchamber)  from  chamber  cross  sectional  area 


D, 


chamber 


=2 


^chamber 


71 


i 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 


Section  2:  Approach  and  Analysis/Analytical  Process 


General  Algorithm-Chamber  Length 


•  Calculate  drop  evaporation  time  (tD)  from  drop  diameter  and  evaporation  constant 


•  Calculate  drop  distance  traveled  (LD)  from  drop  evaporation  time  and  sheet  velocity 
-  Drop  velocity  is  presumed  equal  to  sheet  velocity 


•  Calculate  chamber  length  required  (Xj)  for  impingement  stream  from  chamber  diameter  and 
impingement  half  angle 


COS0 


•  Sum  the  impingement  stream  distance,  impingement  sheet  length,  and  drop  distance  traveled  to 
estimate  required  chamber  length  for  a  given  configuration 


Lr=Xj+rb+LD 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  For  hydrazine,  6  cases  did  not  meet  the 
requirements  of  the  relationships  set  forth  by 
Ryan  et  al.1 

-  Case  2.1.2-2.3.2  and  3.1.2-3.3.2 

•  We  exceeded  upper  boundary 

•  Injection  velocities  were  highest  of 
cases  considered 

-  Uinj  >32  m/s 

-  Smaller  d0  values  relative  to 
Ryan  and  increasing  Pc,  push 
increase  in  Uinj 

-  Ryan’s  highest  Uinj  was  18.5 
m/s 

-  Lowest  N2H4  case  values  was 
15.9  m/s-all  others  exceeded 

-  pL  to  ctl  value  of  N2H4  lower 
than  H20,  however, 
substantially  increased  Uinj 
drove  We  values  out  of  range 


N2H4  @  5  lbf:  We  as  a  function  of 


■  We  lower  =  350 

- We  upper  =  6600 

—♦—Case  1.1. 1-1.3. 1 
—♦—Case  1. 1.2-1. 3. 2 
-■-Case  2.1. 1-2.3. 1 

■  Case  2. 1.2-2. 3. 2 
-•-Case  3.1. 1-3.3. 1 
—•—Case  3. 1.2-3. 3. 2 


Uinj  (m/s) 


•  All  cases  met  Re  requirement 

AFRL 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•For  hydrogen  peroxide,  12  cases  did  not  meet 
the  requirements  of  the  relationships  set 
forth  by  Ryan  et  al.1 

-  Case  1.1.2-1.3.2, 2.1.2-2.3.2, 3.1.1-3.3.1 
and  3.1.2-3.3.2 

•  We  exceeded  upper  boundary 

•  Injection  velocities  were  highest  of 
cases  considered 

-  Uinj>27  m/s 

-  Smaller  d0  relative  to  Ryan 
and  increasing  Pc  push 
increase  in  Uinj 

-  Lowest  H202  case  value  was 
19.1  m/s-greater  than  Ryan’s 
highest  value 

-  Combination  of  increased  pL 
to  ctl  value  of  H202  relative 
H20  &  increased  Uinj  drove 
We  values  out  of  range 


H202  @  5  lbf:  We  as  a  function  of 


■  We  lower  =  350 

■  We  upper  =  6600 
— *~Case  1.1. 1-1.3. 1 

Case  1.1. 2-1. 3. 2 
— B —  Case  2. 1.1-2. 3.1 
— B — Case  2. 1.2-2. 3. 2 
— ©—Case  3.1. 1-3.3. 1 
—©—Case  3. 1.2-3. 3. 2 


•  All  cases  met  Re  requirement 


1 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  Relative  to  Ryan  et  al.1,  the 
values  of  impingement  sheet 
length  to  injector  orifice 
diameter  are  low 

-  Closest  agreement  observed 
for  the  turbulent  cases  of 
empirical  data 

•  These  study  points 
had  the  lowest 
injection  velocities 

•  Larger  values  of  calculated  Uinj 
in  the  study  drive  increased  We 
&  rb  values  relative  to  the  model 
of  Ryan 

-  Case  1.1. 1-1.3. 1  closest  in 
agreement 

-  Increasing  We  of  study  does 
not  overcome  effect  of 
decreased  rb  to  d0  ratio 


rb/d0  versus  We(l-cos0)2/sin30 

Instant  study  results  (20=80°,  do=0.39-0.51  mm)  compared  to  Ryan  et  al.1 

(d  o=0. 64  mm) 


— ♦—  N2H4  cased. 1-1.3. 1 
— ♦—  N2H4  case  1.1. 2-1. 3 .2 
— N2H4  case2.1. 1-2.3. 1 
— N2H4  Case  3.1. 1-3.3. 1 
—  H202  Case  1.1. 1-1. 3.1 
~B—  H202  Case  2.1. 1-2.3. 1 


I  w 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  Relative  to  Ryan  et  al.1,  the  values 
of  evolved  droplet  to  injector 
orifice  diameter  are  low 

-  Appear  to  agree  in  trend 

•  Hydrazine’s  physical  property 
values  are  closer  to  water  (liquid 
Ryan  used)  than  hydrogen 
peroxide 

•  Instant  study  calculated  droplet  to 
orifice  diameter  ratio  appears  to 
be  about  a  factor  of  3  less  than 
those  of  Ryan 

•  The  value  of  $  (3i  mdt~26  for  the 
instant  study,  calculated  with 
physical  properties 

•  Smaller  d0  study  values  resulted 
in  smaller  values  of  dD,  pushing 
study  values  of  their  ratio  below 
those  of  Ryan 


dD/d0  versus  We(l-cos0)2/sin30 
Instant  study  results  compared  to  Figure  8 

of  Ryan  et  al.1 
20=80°,  do=0.39-0.51  mm 


— ♦ — N2H4  Case  1.1. 1-1.3. 1 
— ♦—  N2H4  Case  1.1.2-1.3.2 
— ■— N2H4  Case  2. 1.1 -2.3.1 
— N2H4  Case  3.1. 1-3.3. 1 
-- H202  Case  1.1. 1-1.3. 1 
— S —  H202  Case  2.1. 1-2.3. 1 

0  300  600  900  1200  1500 

We(l-cos6)2/sin36 


1 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


Chamber  dimensions  decrease  as 
the  chamber  pressure  increases 


-  For  a  fixed  thrust  level  this 
is  sensible  since  the  same 
amount  of  high  temperature 

D* 

gas  is  evolved 

600 

•  The  only  path  to 
increased  chamber 

500 

pressure  is  decreased 
volume 

400 

Pc  (psia)  300 

Hydrogen  peroxide  values  are 
slightly  larger  than  those  for 

200 

hydrazine 

100 

-  Lower  specific  impulse 
increases  the  mass  flow  rate 
required 

-  Has  a  higher  decomposition 
temperature 

-  Both  of  these  factors  drive 
increased  values  relative  to 
hydrazine 

0 

D*,  D  as  function  of  Pc  &  y  values  for  N2H4 


&  h2o2 


^ - 


• 


j  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 


•  N2H4  D* 

♦  H202  D* 
— • — N2H4  D 
— ♦ — H202  D 


1  2  3  4  5  6  7 

D*,  D  (mm) 


9  10  11 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  Increasing  chamber  pressure  shows  a 
decreasing  chamber  length  result 

-  Consistent  with  diameter  and 
throat  calculations 


•  Decreasing  1st  stage  mass  flow 
(increasing  second  stage  mass  flow 
rate)  was  not  a  driver 

-  Holding  CD  constant  held  the 
injection  velocity  constant, 
adjusted  the  orifice  dimension 
(increased  with  increasing  mass 
flow  to  the  second  stage) 


N2H4  Chamber  Dimensions 


B 

10.00 

9.00 

S, 

8.00 

5m 

OJ 

« 

7.00 

s 

5 

6.00 

« 

o 

5.00 

u 

J3 

H 

4.00 

*3 

5m 

3.00 

QJ 

pQ 

s 

a 

2.00 

u 

1.00 

0.00 

4.5  9.5  14.5 

Chamber  Length  (mm) 


— N2H4  Case  1.1. 1-1.3. ID 
chamber 

— N2H4  Case  1.1.2-1.3.2D 
chamber 

— ■—  N2H4  Case  2.1. 1-2.3. ID 
chamber 

— •—  N2H4  Case  3.1. 1-3.3. 1  D 
chamber 

— ♦—  N2H4  Case  1.1. 1-1.3. ID 
throat 

— ♦—  N2H4  Case  1.1.2-1.3.2D 
throat 

— *—  N2H4  Case  2.1. 1-2.3. ID 
throat 

— N2H4  Case  3.1. 1-3.3. 1  D 
throat 


•  Increasing  hardness  ratio  decreases 
chamber  length 
-  Increases  injection  velocity 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  Similar  trends  as  with  hydrazine 


•  Increasing  chamber  pressure  shows  a 
decreasing  chamber  length  result 

-  Consistent  with  diameter  and 
throat  calculations 


•  Increasing  hardness  ratio  decreases 
chamber  length 

-  Increases  injection  velocity 

•  Resulting  overall  chamber  length  is 
greater  than  that  of  hydrazine 

—  Greater  relative  injection  velocities 
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a  10.00 

E 

U 

QJ 

ts 

S 

cs 


3.00 


cs 

o 

5- 

pC 

H 

QJ 

pD 

s 


6.00 


4.00 


2.00 


0.00 


H202  Chamber  Dimensions 


m 


4.5  9.5  14.5 

Chamber  Length  (mm) 


19.5 


H202  Case  1.1.1- 

1.3.1  D  chamber 

H202  Case  2.1.1- 

2.3.1  D  chamber 

H202  Case  1.1.1- 

1.3.1  D  throat 

-S — H202  Case  2.1.1- 

2.3.1  D  throat 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  3:  Results  and  Discussion 


•  Contrasting  hydrogen  peroxide  to 
hydrazine,  an  increased  chamber 
volume  is  required  for  an  equivalent 
operating  condition 

-  Most  likely  driven  by  the  decreased 
energy  content  which  drives 
increased  propellant  mass  flow 
rates 


Chamber  Dimensions,  H202  &  N2H4 


Chamber  Length  (mm) 


— o — 
— ■ — 
— ^ — 
— H-- 


N2H4  Case  1. 
chamber 
N2H4  Case  1. 
chamber 
N2H4  Case  2. 
chamber 
N2H4  Case  3. 
chamber 
N2H4  Case  1. 
throat 

N2H4  Case  1. 
throat 

N2H4  Case  2. 
throat 

N2H4  Case  3. 
throat 

H202  Case  1. 
chamber 
H202  Case  2. 
chamber 
H202  Case  1. 
throat 

H202  Case  2. 
throat 


1.1- 1. 3.1  D 

1.2- 1. 3. 2  D 

1. 1- 2.3. ID 

1. 1- 3.3. ID 

1. 1- 1.3. ID 

1.2- 1. 3. 2  D 

1. 1- 2.3. ID 

1. 1- 3.3. ID 

1. 1- 1.3. ID 

1. 1- 2.3. ID 

1. 1- 1.3. ID 

1. 1- 2.3. ID 


I V 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Section  4:  Summary  and  Conclusions 


•  Performed  study  to  develop  first  iteration  chamber  dimension  estimates  based  upon  propellant 
physical  properties 

-  Chamber  dimensions  derived  from  atomization  characteristics  of  propellants  of  interest 

-  Performed  relative  comparison  of  results  to  published  empirical  data 

*  Results  had  small  overlap  with  published  data 

•  Increasing  chamber  pressure  and  hardness  ratio  decrease  the  overall  chamber  dimensions  for  a  given 
thrust  level 

•For  the  propellants  considered,  energy  content  appears  to  have  a  greater  impact  on  the  chamber 
dimensions  than  the  physical  properties 

-  Injection  velocity  appeared  to  be  a  greater  contributor  to  We  number  than  physical  properties 

-  Increasing  We  for  a  given  chamber  pressure  increases  chamber  length  required  to  achieve  liquid 
to  vapor  transition  of  droplets 

•  Relative  comparison  with  resulting  engine  dimensions  from  study  difficult  to  compare  to  state  of  the 
art  engines 

-  Requires  knowledge  of  injection  parameters,  generally  held  proprietary 

-  By  visual  inspection,  the  dimensions  are: 

•  Small  relative  to  state  of  the  art  monopropellant  hydrazine  &  nitrogen  tetra- 
oxide/hydrazine  bipropellant  engines  for  the  same  thrust  level 
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Nomenclature 


Notation 

Definition 

a 

Curve  fit  polynomial  constant 

A  chamber 

Chamber  cross-sectional  area 

A* 

Throat  cross-sectional  area 

A0 

Injector  orifice  cross-sectional  area 

Bq 

Spalding  or  transfer  number,  based  on  heat  transfer  considerations  only 

CPg 

Specific  heat  at  constant  pressure  of  gas  phase 

Cd 

Discharge  coefficient 

CDmax 

Discharge  coefficient  value  attained  at  Re  >  10,000 

CH3OH 

Methonal 

cp 

Specific  heat  at  constant  pressure 

Cv 

Specific  heat  at  constant  volume 

Dchamber 

Chamber  diameter 

D* 

Throat  diameter 

<Id 

Droplet  diameter 

d0 

Injector  orifice  diameter 

djet 

Impinging  liquid  propellant  jet  diameter 

Ft 

Thrust  force 

go 

Gravity 

h 

Liquid  propellant  sheet  thickness 

H 

Hydrogen,  monatomic 

h2 

Hydrogen,  diatomic 

h2o 

Water 

h2o2 

Hydrogen  Peroxide 

hfg 

Latent  heat  of  vaporization 

Isp 

Specific  impulse 

K 

Evaporation  constant 

kp 

Propellant  gas  phase  thermal  conductivity 

kg 

Environment  gas  phase  thermal  conductivity 

km 

Wave  number  for  the  most  unstable  wave 

kinf 

Gas  generator  combustion  product  gas  phase  thermal  conductivity 

Ld 

Distance  traveled  by  droplet  until  fully  evaporated 

Lt 

Total  chamber  distance  required  for  injection,  atomization  and  evaporation  of 
liquid  propellant 

lo 

Injector  tube  length 

M 

Mach  number 

m’ 

Mass  flow  rate  of  liquid  propellant 

MeOH 

Methonal 

m’gg 

Gas  generator  mass  flow  rate  of  liquid  propellant 

m’gg/m’inj 

Ratio  of  mass  flow  rate  of  liquid  propellant  injected  through  the  gas 
generator  to  the  that  injected  downstream 

m’inj 

Injected  mass  flow  rate  of  injected  liquid  propellant 

m’T 

Total  mass  flow  rate  of  liquid  propellants  injected  into  the  thruster,  derived 
from  specific  impulse  and  thrust  requirement. 

MSDS 

Material  Safety  Data  sheet 

MW 

Molecular  weight 
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Impact  of  Liquid  Propellant  Properties  on  Small  Rocket  Thruster 

Dimensions 

Nomenclature 


ATjnf.boil 

The  difference  between  the  free  stream  environment  temperature  and  the 
liquid  propellant  boiling  point 

Uinj,  Uj 

Injected  liquid  propellant  jet  velocity 

Us 

Velocity  of  liquid  in  sheet  resulting  from  impingement  of  injected  liquid  jets 

Vac 

Vacuum 

We 

Weber  number,  ratio  of  inertia  to  surface  tension 

X 

Percentage  of  ammonia  dissociation  in  hydrazine  decomposition  reaction 

xj 

Horizontal  length  component  of  liquid  jet  within  chamber 

Yi 

Molar  mass  fraction, 

Pi,m 

Liquid  sheet  disturbance  maximum  growth  rate 

Y 

Specific  heat  ratio 

Ymix 

Specific  heat  ratio  for  a  particular  molecule 

Pl 

Dynamic  viscosity  coefficient 

Vl 

Kinematic  viscosity  coefficient 

♦ 

Angle  relative  to  centerline  of  the  liquid  sheet  formed  from  impingement  of 
liquid  propellant  jets 

PL 

Liquid  propellant  density 

Pc 

Chamber  gas  phase  density 

oL 

Surface  tension 

0 

Impingement  half  angle 

AFRL^ 


N 

Nitrogen,  monatomic 

n2 

Nitrogen,  diatomic 

n2h4 

Hydrazine 

ni 

Mol  fraction  for  a  particular  species 

nx 

Total  mols  of  all  species  in  a  reaction 

0 

Oxygen,  monatomic 

o2 

Oxygen,  diatomic 

Pc 

Chamber  pressure 

Pf 

Feed  pressure 

Ptfc 

Ratio  of  feed  to  chamber  pressure,  other  wise  referred  to  as  “hardness”  ratio 

rb 

Radial  distance  from  impingement  point  where  liquid  sheet  break  up  into 
ligaments  occurs 

R,RU 

Universal  gas  constant 

Rspec 

Universal  gas  constant  divided  by  the  MW  of  the  gas  of  interest 

Re 

Reynolds  number,  ratio  of  inertial  forces  to  viscous  forces 

s 

Ratio  of  chamber  gas  density  to  liquid  propellant  density 

SL 

Sea  Level 

T 

Temperature 

Tbar 

Average  temperature  of  the  free  stream  environment  and  the  liquid  propellant 
boiling  point 

Tboii 

Boiling  point  temperature  of  liquid  propellants 

Tc 

Chamber  (combustion  product)  temperature 

Tinf 

Free  stream  environment  temperature 
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